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Abstract 

Sequences of canonical conservation laws and generalized symmetries for the lattice Boussinesq and the lat- 
tice modified Boussinesq systems are successively derived. The interpretation of these symmetries as differential- 
difference equations leads to corresponding hierarchies of such equations for which conservation laws and 
Lax pairs are constructed. Finally, using the continuous symmetry reduction approach, an integrable, multidi- 
mensionally consistent system of partial differential equations is derived in relation with the lattice modified 
Boussinesq system. 

1 Introduction 

The existence of infinite hierarchies of symmetries and conservation laws is among the key characteristics of most 
of the well known integrable systems, (7] Q21 1201. Recently, in (8), this characteristic of integrability was also 
established for a class of two-dimensional integrable partial difference equations, namely the equations in the 
Adler-Bobenko-Suris (ABS) list, cf. UJ, of canonical scalar affine-linear lattice equations. Moreover, in |[T8l . it 
was shown that equations in this class possess infinite hierarchies of symmetries and conservation laws, and these 
hierarchies were constructed recursively. Their derivation is similar to the approach used in 0, and it is based on 
the use of the first canonical conservation law and the first master symmetry. In the present paper, we apply the 
same ideas to the case of the lattice Boussinesq and lattice modified Boussinesq equations. 

By the lattice Boussinesq equation^, we mean the following partial difference equation on the two-dimensional 
lattice 

a 3 - f3 3 a 3 - f3 3 

7T~ — 7r~ — Mn,m+l M n+l,m+2 + «n+l,m' i n+2,m+l 

OL — p + U n+ i^ m+ i — U n+ 2, m a — p + U n ^ m+ 2 — U n+ i^ m+ i 

+U n +2,m+2 (a — /3 + U n+ i jm+ 2 — U n +2,m+l) + Un,m (" ~ P + Un,m+1 ~ Un+l,m) 
= (2a + P) (Un+l,m + Un+l,m+z) ~ ( a + 2 /3) («n,m+l + U n+2 ,m+l) , (1) 



and we refer to 



OC 2 U n+ i^ m+ i — (3 2 U n ^ m+ 2\ n n+ i j?n+ 2 / a 2 U n+ 2,m — P 2 U n +l,m+l \ ""n+2,m+l 



a^n,m+2 — pUn+l,m+l J ^n,m+l \ (™n+l,m+l — Pu n +2,m J ^n+l,m 

a \ (2) 



, u n+l,m U n +2,m+2/ \Un,m+l ^n+2,m+2 , 

'An alternative version of a lattice Boussinesq equation was given in |2| which could be viewed as a dimensional reduction of the 
modified lattice KP equation. However, it seems that this equation corresponds to the special case p = xi, where \ is a primitive root of 



unity, of equation ([TJ. 
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as the lattice modified Boussinesq equation. These equations were derived in Ifl2ll in the context of the discrete 
Gel'fand-Dikii hierarchy, and, under appropriate continuum limits, they reduce to the continuous potential Boussi- 
nesq equation and the potential modified Boussinesq equation, respectively. These 9-point scalar equations can 
also be given as quadrilateral systems for three or two fields, respectively, and these forms will be used in the 
subsequent analysis. 

In particular, our analysis will be focused on the derivation of hierarchies of symmetries and canonical conser- 
vation laws for these systems. These hierarchies are derived in a systematic way starting with the first generalized 
symmetry and the first canonical conservation law as seeds and making successive use of an extended generalized 
symmetry playing the role of a master symmetry. This procedure was used in [ 1 8 ] in the derivation of correspond- 
ing hierarchies for all the ABS equations. 

Moreover, for the system related to the lattice modified Boussinesq equation, a corresponding system of partial 
differential equations (PDEs) is derived similarly to the case of the ABS equations Ifl6l . These systems are also 
referred to as generating PDEs and were introduced in IPTT1 in connection with the lattice potential KdV equation. 
Subsequently, similar systems derived for the lattice Boussinesq equation [ 14] and for all the ABS equations ifloTl . 
These systems inherit some properties of their lattice counterparts, more precisely they admit an auto-Backlund 
transformation and a Lax pair, both of which actually follow from their multidimensional consistency. For the 
multidimensionally consistent generating PDE related to the lattice modified Boussinesq equation, we also present 
an auto-Backlund transformation and a Lax pair. 

The paper is organized as follows. The next section deals with the symmetries and conservation laws of the 
lattice Boussinesq system, while Section[3]contains the corresponding results for the modified Boussinesq system. 
In Section |4] the generating PDE related to modified Boussinesq system is presented along with an auto-Backlund 
transformation as well as a Lax pair for this system. The concluding section discusses some related points of 
current and future research. 



2 Lattice Boussinesq equation 

In this section we study the generalized symmetries and conservation laws of the lattice Boussinesq equation. 
This equation was derived in Ifl2l . studied in [ 17] while the corresponding generating PDE was presented in Ifl4l . 
Here we derive formulas for the recursive construction of symmetries and conservation laws using as seeds a 
pair of symmetries and a canonical conservation law. The interpretation of these symmetries as a hierarchy of 
differential-difference equations and our considerations lead naturally to a hierarchy of conservation laws and Lax 
pairs for these hierarchies as well. 

Let us begin by writing the lattice Boussinesq equation as a system on a quadrilateral. Specifically, this equa- 
tion can be written as a system involving the values of three fields u n m := (u njm , i> n ,m, ^n,m) on an elementary 
quadrilateral of the lattice which has the following form. 

^n+X, m V"n,m'U"n+l,m ~\~ Vn,m — > (3a) 
W n ,m+l "n,m%,m+l Vn,m — > (3b) 

a — (3 

'Un,m'U'n+l,m+l ^n+l,m+l ^n,m — . (3c) 

A Lax pair, gauge equivalent to the one derived in lfl2l . cf. 0/71 [14]], is given by the following linear system for 

0) V n ,m , (4a) 
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where 

L(u niJn ,u„ + i im ;a) := — _^ 1/3 | -«n+l,m 1 ) . (4b) 

The compatibility condition of the above linear system holds if and only if system © holds. 

One useful property of system (0) is its covariance, i.e. its invariance under the mutual interchanges of shifted 
values of u, v and w accompanied with corresponding interchange of the parameters a and (3. This property, which 
is obvious for the Lax pair ©, allows us to present symmetries and conservation laws only for one lattice direction, 
specifically for the n direction, while the results for the m direction can follow by employing the interchanges 

(^71+1,771+^ ) ^Vi+i,m+j) Wra+i,m+j) HI, Ct, /3) \ — (Un+j^n-j-j-j, V n -\-j )m + j , U) n _|_j iJn -|_j , 771, 77, /3, tt). 

Theorem 1 J/ze lattice Boussinesq system (0 admits a ten-dimensional algebra of Lie point symmetries generated 
by the vector fields 

So du m + 

*V n m 

9v n ,m + 2w n,mdw n , m + 3a<9 Q + 3f3dfs , 
t^o = d Un m + u n>rn d Vn m + u n ^ m d Wn m , 
To = 9i) n m — d Wn m , 

£7, — -v k ( n + m )v 8 — -yHn+m+i) o _ fc(n+m+l) a 

rr, _ ~.k(n+m) f) I ^.fcfn+m+l),, i ,.fc(n+m+2)„. a 

^fc — A ^Tin.m ' A u n,m u v„, m i A ii n,m'Jw n>m > 

T, — - v/ fc(7l+77l) a _ , / fc(77+777+2) 1 f) 

A k — A U V ntm A U W n ,m 1 

E = d a + dp, 

where k = 1, 2 ant/ % jj a primitive root of unity x 2 + X + 1 = 0. 

Moreover, it admits a hierarchy of generalized symmetries Vi and an extended generalized symmetry G\ in the 
n direction, which have the following form 

V* = ^ ( —) d Un , m + ( + ( ^=m) flL , 7 = 1,2,-.., (6) 



and 

Gi = d Un m + d Vnm + d w + d a . (7) 

' 77,771 ' 71,777 ' 77,771 

Function r n>m and operator ffl are given by 

Tn,m ■ "^77+1,777^77—1,777 ^77+1, 777 ^77 — 1,777 (8a) 

and 

OO oo oo 

J O , 3 77+1+^,777 o J ^77 — 1+^,777 



*== E E =V^Wi,~ + I! J -f=^0 Wn+ , m + d a , (8b) 

'77+1,777 • '77+7,777 • ' 77 + 7,777 

J= — OO 1 ■" J = — OO "" J = — OO J ' 

Moreover, system <\3j admits a hierarchy of canonical conservation laws A m pn,m = &-n&n,m, the densities 
Pn)n and the fluxes o~H%, of which are defined recursively as 

Pn)n = ( l0 S (^77+1,777^77-1,777 ~ ^77+1,777 ~ W n -l >m ) ^ , (9a) 

a i k m = ( lo g (Un,m+lU n -l,m ~ V n ,m+1 ~ W n ^ m ) 2 ^ , (9b) 
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for k = 0, 1, • • •. Operator ffl* is the "negative part" of&, i.e. 

®* ■= E -^^ +Jim + g ^±^^„ + , m + £ ^^^ + , m + a Q . (io) 

1 n+j,m . 1 n+j,m ■ ' n+?',m 

j=—oc •" j=~oo ■" j=—oo 

Proof It can be verified directly that the vector fields given in (|5]> are symmetry generators for system ([3]). They 
form a ten dimensional solvable Lie algebra with the non-vanishing commutators given by 

[S Q ,Ui] = -Ui, [S k ,Ui] = -U k+ i mod 3, 
[So, Ti] = -2 T t , [S k , T e ) = X 2k T k+i mod 3 , 
[U k ,Ut] = (X e - X k )T k +l mod 3, 
[S ,E] = -3E, 

where k= 1,2 and I = 0, 1,2. 

For the hierarchy ©, we observe that G\ is a master symmetry of V\. Indeed, their commutator is non-zero, 
resulting to V%, and [Vi, [G±, Vi]] = 0. This implies that we can construct recursively the whole hierarchy © of 
commuting symmetries by setting Vi + \ = [G\, Vi]. Working out these commutators, relations © easily follow. 

Having constructed the hierarchy of generalized symmetries, one may derive hierarchies of canonical conser- 
vation laws starting with the first canonical conservation law A rn p^m = A n o"i°m. The density pn)n and the flux 
Un% have the form 

Pn°L = lo g r n,m and a n)n = lo S («n,m+l«n-l,m ~ ^n,m+l ~ W n -l,m) 2 , (Ha) 

respectively. The form of the conservation law suggests the introduction of the potential (j> (°) through the relations 

A n 0(°)=p(%, A n ^=a£l. (lib) 

If we consider the system constituted by © and (fTTb . then its lower order generalized symmetries follow by 
extending the action of V\ and G\ in the direction of the potential (j> (°) . This means that they have the form 

t/ _ 1 a i u n+l,m a U n -\ %m i(f,(°)« 

V l - Vun.m H ° F «n,m H + *1 °U (0) > 

' n,m ' n,m 1 n,m 

and 

Gi = 3 u „, m + —d v + '— d Wn<m + (n*i + w )fyo» + d Q , 

' n,m ' n,m 1 n,m 

respectively. The condition that these are symmetry generators for system (fTTT> implies that 

*(0) _ xr (AO) \ = (J_f) | "n+l,m a , "n-l,m » V (0) 

v 7 \'n,m I n,m 1 n,m / 

and W is the potential corresponding to the next conservation law. In particular, the density pn) m and the flux 
&n]m of this conservation law are given by 

p£l = & (pi l) - ($S 0) ) , 4]l = ®* {4%) d2) 



where is the "negative part" of St given in (fTOl . Moreover, taking into account the particular form of ^\ 
density pn)n may also be written as pi]m = (pn% 
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If we include potential <j> W i n our considerations and extend symmetries V\ and G\ to this direction, we will 
introduce a potential (f> ^ which corresponds to the next canonical conservation law. We can follow the same 
procedure with this new potential which will result to the introduction of another potential and corresponding 
conservation law. This procedure can be repeated infinitely many times, leading to a hierarchy of conservation 
laws. In particular, symmetry generators will have the form 

t r 1 o . u n+l,m n U n —i m \ „ 

Vi = d Un ^ m + d Vnim + d Wn>m + 2^ *i , 

Tn,m Tn,m Tn,m , 

G l - "lin.m H 0, «n,m ^ ~ °W,m + / A n$ l + V J°V>(*0 + "<* ! 

<n,m 'n,rn in,m , 

where 

^1 - ^1 I /V-l,m ) = I ~ d un, m + — (J Vn m + — d Wn m p n _ 1>m . (13) 

^ / \'n,m ' n,m ' n,m J 

(k) a (k) 

Using the last relation, the corresponding densities and fluxes of the conservation laws A m pn,m = A n <Tn,m can 
be written as in d9ab . 

Moreover, it is straightforward to derive the action of V,, i > 1, on the potentials using the definition of 
these symmetries and the above extended forms of V% and G\. This derivation results to 

which is valid for i = 1,2, • • •, and 

i=o V J / 

Using relations ([T3l and (|9al ), the above relation can be written also as 

= E ( ' ~ 1 ) (-l)^ 3 - 1 ® 3 o Fx o J?"* o^-H-i" 1 (logr n 2 m ) , (14b) 
j'=o 

which hold for i = 1, 2, ■ • ■ and fc = 0, 1, 2, ■ ■ •. □ 



The interpretation of symmetries as differential-difference equations commuting with the lattice system and 
the previous analysis lead to hierarchies of differential-difference equations and corresponding conservation laws. 
Since they involve shifts only in one lattice direction, in what follows we drop the second index m which does not 
vary. In this setting, we have the following 

Corollary 1 The hierarchy of differential-difference equations 

dU \r n J dti V r n J ' dti V r n 

where 

r n := u n+1 u n -i - v n+1 - w„-i (16a) 



1,2,- 



(15) 
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and 

CO CO CO 

* := £ + E + E + ft. . (1*) 

' n+? . ' n+ j . ' n+i 

j=— co J j=— co J j=— co J 

admits a master symmetry 

dun _ _n_ dvn _ _ nun-i n ~ 

oa r n oa r n oa r n 

Moreover, systems f lTTl ) admit a hierarchy of canonical conservation laws, which have the following form 

|:(^ fe (logr n 2 )) = A n (£( ^ ) (-l^-^oro^o^^-^logr^j , 

w/f/z = 0, 1, 2, • • •, i = 1, 2, • • •. In the above relation, is the "negative part" of ' ffl given in M6b\ and 

V — — d I Un+l d I Un - l d 

riy* <Y* 
n 1 n 1 n 

Finally, a Lax pair for systems f li 51) and 47 71 ) is given by 

^ n+ i = L(u n ,u n+1 ;a)fy n , (18a) 



da 

Matrix L is defined in H4b\l and traceless matrix AfH is given by 
Af(u n ,u n+ i,u n _i;a) : = 



M(u n ,u n+ i,u n _i;a)) ^ n , i = l,2, (18b) 
nM(u n ,u„ + i,u n _i;Q!)^ n . (18c) 



%+l-2iUn-l-«n-lM?i + l 1 

u n— 1 1 



11/ 3 

-Mn+1%-1 3 — Un+1 | • UoOJ 

„, „„ „, „. m„-l-2tl„ + l-M ra -llt, 1 + l 

—Vn+lWn-1 U n -lV n+ i , 



a — A TV, 



Here we have constructed infinite sequences of symmetries and conservation laws for system ® which is 
related to the lattice Boussinesq equation £[]). In the next section we are going to establish similar results for the 
lattice modified Boussinesq equation © written as a system for two fields. 



3 Lattice modified Boussinesq equation 

The modified Boussinesq equation was derived in Ifl2l and studied in ifTTl . while its Lie point and lower order 
generalized symmetries were presented in [15]. In this section we present hierarchies of symmetries and canonical 
conservation laws for this equation, as well as corresponding structures for differential-difference equations. We 
follow the derivation of the previous section and, thus, we omit most of the details here and present only the 
necessary results. 

2 The continuous part d 1 8cb of the Lax pair for the master symmetry dl7t is gauge equivalent to the one derived in 1 14 1. 
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The lattice modified Boussinesq equation can be written as a system for two fields u n?n := (ii n , m , f n ,m) 
defined on an elementary quadrilateral of the lattice in the following form 

— /3u n+ i m 

U n +l,m+l = V ntm , (19a) 

<XVn+l,m — PV n ,m+l 

Vn,m Q ; 'f , n+l,m'^n,m+l /^7i,77i+l'^7i+l,ra , lnl . 

V n +l,m+l = — ! ! ' —3 — 1 — , (19b) 

^n,m Q^n+l,m P^V^m+l 

cf. |[T0l where a similar 2-field system for the lattice modified Boussinesq equation was first given. This system 

is the necessary and sufficient condition for the consistency of the following linear system, which, in this sense, 
constitutes a Lax pair for system ( fl9l) . 



,m+l — ^ J \}^-n,mi U n ,m+1 ;/3)*n,m, (20a) 

where 



L(u ri:m , u n+ i jm ; a) 



(a 3 - A 3 )V3 



^ Qttn+l,m^n,m 



V 

The point symmetries of this system, which were given in [15 ], 



S[X) = X n+m Un,md Un>m ~ X n+m+1 Vn,md Vn , m , (21a) 



where \ is a primitive root of unit, \ 2 + X + 1 = 0. Additionally, system (fT9l admits an extended point symmetry, 
namely 

£ = ad a + . (21b) 

Theorem 2 77ie modified Boussinesq system < \19i admits a hierarchy of generalized symmetries V{ and an extended 
generalized symmetry G\ which have the form 

V t = & i - 1 (R 1 )a Un<m + &' i - 1 (Pi)9v n , m , i = 1,2,-.., (22a) 

and 

Gi=n Rid Un m + n P\d Vn m - ad a , (22b) 

respectively. In the above relations, 

Hi := u n , m , (22c) 

Un-\-l,m'^n,m ~r "n,m"n-l,m t %-l,m^n+l,m 

A := 1 —- 1 h w n , m , (22d) 

^n+l,m,Vn,m ~T~ Wn,m"n-l,m T W-n— 1,771^71+1,771 

ant/ 

oo oo 

^ := E i^'(^iR n+J , m + E j^i(Pl)dv n+j>m -ad a . (22e) 

j = — OQ j = ~OD 
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(k) a (k) (k) 

Moreover, it admits a hierarchy of canonical conservation laws A m pn,m = An&n,m, the densities pn,m and 
(k) 

the fluxes o~n,m of which are given by 

p W =(^ + 9] k ( log ( ^±l^IbB + + Un-l,mV n+l , m \ 2 \ ; (23a) 



(fc) _ cpk ( i ( a2u n,mVn-l,m + a/3Un,mVn,m+l + P 2 Un,m+lVn,m X 2 



(23b) 



for k = 0, 1, 2, ■ • ■. f/ze above relations, operators £F a«<i are given by 
and 

-l -l 

^ := E i^(^i)^n + „ m + E i^( P l)^n + ,-, m -«^ S (24) 

j=—oo j=—oo 

respectively. 

Proof Using the fact that Gi, given in (I22bl ). is a master symmetry for V\, hierarchy (I22ab can be constructed 
recursively by setting V^ + i = [G±, Vj\ for i = 1, 2, • • •. 

The hierarchy of conservation laws with densities and fluxes given in (I23ab can be constructed using the same 
method and analysis we used in the previous section. That is starting with 

Pl/m = log ( J , 

\ Uj n,m u n,m / 

(0) 1 l a u n,mVn—l,m ^71.771+1^71,771 \ 

and employing the first symmetry V% and its master symmetry G\. In the course of the construction, we extend 
symmetries V\ and Vi,i > 1, in the direction of the corresponding potentials and these extensions are given by the 
functions 

*?° = - 2u n , m )d Un , m + (P 1 - v n>m )d Vn>m ) pW ltm , k = 0, 1, • • • , (25a) 

and 

$ (fc) = / i - 1 \ ^y+j-i^j ^(i+fe-j-i)^ ^ fe = 0, 1, • • • , (25b) 

respectively. □ 
Remark It is worth mentioning that the first generalized symmetry V\ is actually a linear combination of two point 
symmetries and a generalized one, specifically 

rj o I Q I ^71,771 "^71+1,771 ^71, 771 q ^71,771 ^Tl— 1,771^71,771 o 

"l ~~ ~ u n,mC'u n , m ~r V n>rn O Vn m + a I "Wn.tra "> "t> n , 

\ 'n,m Tn,m 

Our choice for the form of Vi is motivated from the corresponding form of Gi and does not affect our analysis. 
The interpretation of symmetries as differential-difference equations leads to 
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Corollary 2 The hierarchy of differential-difference equations 



dUr, 

dti 



admits a master symmetry 
In the above relations, 



a- 



dun 
da 



dv r , 



-nR\. a 



^-\Pi), * = 1,2,---, 



R, :-- 



Pi :-- 



da 
3u n u n+ iv n 



-nPi . 



Un+lVn + U n V n -l + U n ^iV n+ i 

-3u n v n _iv n 



U r , 



+ V ri 



and 



0> .._ 



Un+lVn + U n V n -l + Un-iV n +l 
oo oo 

E jy£(Ri)d Un+j + E j^i(Pi)dv n+j -cxd 



(26) 
(27) 

(28a) 
(28b) 

(28c) 



j=-oo 



j=-oo 



Moreover, it admits a hierarchy of canonical conservation laws d^pn 1 — A n $-"' the densities of which are 
given by 



p (k) 

rn 



+ 



log 



Un+lVn + UnVn-l + U n -lV n +l 



u n v n 



k = 0,1,2, 



and the corresponding fluxes are defined recursively by 



(*) 



(i?! - 2u n )d Un + (Pi - UnJfl^jpS , 
i-1 



i - 1 
3 



? > 1, 



/or fe = 0, 1, 2, ■ ■ ■. In the above relations, 

Si := 2u n+ id Un+1 + v n+ id Vn+1 
and 3?* is the negative part of & defined in \28c\ . i.e. 



■■= E 3^i{Ri)9 Un+j + E j-^cw 



n+j 



ad n 



j=—oo j=— oo 

A Lax pair for the differential-difference hierarchy \26\ and its master symmetry d27D is given by 

# n+ i = L(u n , VLn+i; a)* n , 

^ = (^^M^Un+i.Un-ija)) * n , 

n , a 3 \ T 
- M(u n ,u n+ i,u n _i;a) + — ^1 ) W n , 



3a a 
where matrix L is defined in H20bi and 

M(u n ,u n+ i,u n _i;a) := 



3 - A 3 



(29a) 
(29b) 

(29c) 



-3a 3 



(a 3 — A 3 )7v, 



XUn+lVn-l 



Un+lVn 



\u n -iv n 
a 

A u n -iu n +iv n 



(30) 



i A Un%- lVn+1 \u n V n +l 
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Thus we have systematically constructed the generalized symmetries and the canonical conservation laws for 
both systems (O and ( fl9l ) which are related to lattice potential and modified Boussinesq equations, respectively. 
The existence of such hierarchies establishes the integrability of these systems, the multidimensional consistency 
of which was established in [17]. In the next section we employ both of master symmetries of system ( fT9l ) in the 
derivation of an integrable and multidimensionally consistent system of PDEs. 



4 Continuous symmetric reduction of the lattice modified Boussinesq system 

The concept of generating PDE was introduced in [ 1 1 ] and a systematic method for the derivation of such systems 
was presented in lfl6l based on the symmetries of the underlying discrete system. Here we apply this method to 
the case of the lattice modified Boussinesq system ( fT9T ) using both of its master symmetries. That is, we look for 
solutions of system (fl9l > which remain invariant under the action of both symmetries, and which lead to a system 
of PDEs rather than a system of ordinary difference equations H. 
More precisely, introducing the notation 

U . — U n m , U\ . — U n +i j?Tl , U2 ■ — "Wn,m+li U—\ . — U n —\^ m , U—2 ■ — U nm — \ , 

v := v n>m , vi := v n +\ >m , v 2 := v n , m+ i, w_i := u n -l,m) V-2 ■= v n ^ m -\ , 

we impose that u, v and their shifts satisfy, additionally to system (fT9l ). the "similarity constraints" 

du ( 3 v u\ \ 
a— \- nu \ 1 = , (31a) 

da \uv-i + u\v + u-\v\ J 

dv ( _ ?>uv-\ \ 
a— + nv 1 =0, (31b) 

oa \ uv-i + u\v + U-\V\ J 

(3— + mu 2 - 1 = 0, (31c) 

OP V™~2 + U 2 V + U_ 2 W2 / 

p — + mv 1 j- =0, (31d) 

Op \ UV-2 + U 2 V + U- 2 V2 J 

which in ifTTl (in a different notation) were presented as a system of differential-difference equations compatible 
with the original lattice system (fT9l ). 

From the above system of differential-difference equations accompanied with the lattice modified Boussinesq 
system (fl9l >. a system of partial differential equations for (u, m, 112) := {u, v, u\, v\, u?, v 2 ) can be derived in a 
systematic way as it was done for the ABS equations in Ifl6l . This system has the following form 

du\ a(au2 — /3ui)(auiV2 — (3u2V\) du f3(aii2 — (3ui)(/3v 2 — av\) dv 

~dp~ = (a 3 - ^)uu 2 v 2 ~dP (a 3 - ^)vv 2 ~d~P 

a(u 2 v 1 (au 2 - Put) + v 2 (Pu\ - au\)) 
—m — ^ -pr- , (j^-a) 

(ce - P 6 )U 2 V2 

dv\ /3{/3v 2 — av\){au\v 2 — /3u 2 v\) du (av\ — j3v 2 ){f3 2 v\ — a 2 V2) dv 
~dp~ = (a 3 - /3 3 )uu 2 V2 ~dfi (a 3 - ^)vv 2 ~d~P 

_ a{uiV2{fiv2 - avi) + u 2 (av\ - /3vf)) 
[a 6 - p A )u 2 V2 



3 We point out that in this respect such symmetry reductions are different from reductions by Lie point symmetries, in that the latter 
usually reduce the space of solutions to a finite-dimensional one, whereas with the reductions we impose here the solution space remains 
infinite-dimensional. 
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du 2 P(/3ui — au 2 ){Pu 2 v\ — au\V2) du a{j3u\ — au2){av\ — j3v2) dv 

da (P 3 — a 3 )uu\Vi da (/3 3 — a 3 )vv i da 

P(uiv 2 {fiui - au 2 ) + vi{au\ - Pup) 

-n 7^3 or , (32c) 

(p J — a 6 )u\Vi 

dv2 a(av\ — Pv2){Pu2V\ — au\V2) du (Pv 2 — av\)(a 2 V2 — P 2 v\) dv 

da (P 3 — a 3 )uu\Vi da (P 3 — a 3 )vv\ da 

f3(u2Vi{avi - (3v 2 ) + uiifivf - gyp) 

d 2 u l_f i+ apr Ul V2 



dad/3 u \ a 3 — P 3 \ U2V\ 



u 2 v\ \ \ du du 
u\V2 J J da dp 



aft av\ — j3v2 ( U2 du dv u\ du dv 



a 3 — (3 3 v \u\V2 da dp u 2 vx df3 da 

a ( ( U2V1 „uiv 2 n u 2 \ du uu\(j3v2 — av\) dv 

~ m_ 3 m a +^ ~ ^~ 

a a _ jj-i U \V2 U2V1 u\ J da U2W1 da 

ft ( ( n u\v 2 U2V1 u\\ du uu 2 (av\ — (3v 2 ) dv 
—n—7z ^ p h a h a 



P 3 — a 3 \\ U2V1 u\V2 u 2 J dp U1W2 dp 

nm (Pv 2 + avi)(u 2 vi — u\v 2 ) 
1 3 — P 3 U1U2V1V2 
d 2 v 1 {a 2 v 2 + P 2 vi)(avi — Pv 2 ) dv dv 



nm (Pv 2 + av^ulvi - ufug) 
a A — p^ uiu 2 v\v 2 



dadp a 3 — P 3 vv\v 2 da dp 

aP au\v 2 — Pu 2 v\ f 1 dv du 1 dv du 



a 3 — P 3 u \u 2 v 2 dadp u\V\dpda 

a ( ( u\ v 2 n vi\ dv v(auiv 2 — Pu 2 v\) du 
-m— 5 -7T a h a h p 



a 3 — P 3 \\ u 2 v\ v 2 J da uu\V\ da 

P ( ( n u i n v i v 2 \ dv v(Pu 2 vi - auiv 2 ) du 
-n—, 3- p hp Va 1 



P 3 — a 3 \ \ u\ v 2 v\ J dp uu 2 v 2 dp 

nm v{auiv 2 + Pu 2 vi)(u 2 v 2 - uivx) 
a A — p^ u\u 2 viv 2 

In the following we will denote this system by i? n m (u, ui, U2; a, P) or, simply, by B n>Tn . 

The integrability properties of this system are inherited to it by its lattice counterpart and are the following 
ones : it is multidimensionally consistent in the sense of lfT9ll , it admits an auto-Backlund transformation and a 
Lax pair. These properties are formulated in the following three theorems which can be proven by straightforward 
calculations. 

Theorem 3 The three copies 

Bn^njiu, Ui, Uj;ai,aj), B nj , nk {u, Uj, u k ; aj,a k ), B nk , ni (u, u k , Uj; a k , a*) 

of the above system are consistent with each other. That is, the three different ways to evaluate d ai d aj d ak u and 
the two different ways to evaluate d ai d a . u^, i ^ j ^ k 7^ i, lead to same results, respectively. 
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X(u 2 v(Xu 2 - f3u) + v 2 (fiul - Xu 2 )) 
m 7-3 ^7 , (33c) 

(X d - P 6 )u 2 V 2 



Theorem 4 Let us denote by B^(u, ui, U2; u, Ui, u 2 ) the algebraic-differential system 

du X(Xui — au)(Xuv\ — au\v) du a(Xu\ — au)(av\ — Xv) dv 

da (A 3 — a 3 )uu\V\ da (A 3 — a 3 )vv\ da 

X(u±v(Xui — au) + v\(au\ — Xu 2 )) 
(A 3 — a 3 )u\Vi 

dv a(avi — Xv)(Xuv\ — au\v) du (Xv — av\)(a 2 v — X 2 v\) dv 

da (A 3 — a 3 )uu\Vi da (A 3 — a 3 )vv\ da 

X(uv\(avi — Xv) + u\(Xv 2 — av 2 )) 
(A 3 — a 3 )uivi 

du X(Xu2 — [3u)(Xuv2 — /3u2v) du (3(Xu2 — (3u)(/3v2 — Xv) dv 
dp = (A 3 - (3 3 )uu 2 V2 d/3 (A 3 - p 3 )vv 2 dp 

dv f3((3v2 - Xv)(Xuv 2 - l3u 2 v) du (Xv - /3ti 2 )(/3 2 -u - X 2 v 2 ) dv 

dp ~ (A 3 - P 3 )uu 2 V2 dp (A 3 - P 3 )vv 2 ~d~P 

X(uv 2 (Pv2 - Xv) + u 2 (Xv 2 - Pv 2 )) 

-m —3 -3- , (33d) 

(A d - P 6 )U 2 V2 

au — Xu\ v au\v — Xuv\ p-u — Xu2 v PU2V — X11V2 

Ul = V — , Vi = -r— , u 2 = v- — , v 2 = -5 r-r— • (33e) 

av\ — Xv u av\ — Xv pv2 — Xv u PV2 — Xv 

System M\(u, u\, u 2 ; u, ui, U2) maps any solution (u, ui, u 2 ) of system B n ^ m to a new solution (u, ui, u 2 ) of the 
same system. In other words, it constitutes an auto-Bdcklund transformation of system A3 '2D . 

Moreover, starting with any solution (u, m, u 2 ) of system B n ^ m and using the auto-Bdcklund transformations 
M\ and to construct two new solutions (u, iii , u 2 ) and (u, ui , u 2 ), respectively, then a third solution (u, iii , u 2 ) 
can be constructed algebraically according to Bianchi commuting diagram, cf. Figure\J\ and this solution is given 
by 

Xu — p,u - v Xuv — uuv 

u = v— -, v = — — , 

Xv — uv u Xv — jiv 

Xii\ — uui t v 1 Xu-ivi — uuiv\ 

ui = vx — — , Vl = — , 

Xvi — uvi u\ Xv\ — uv\ 

z XU2 - UU 2 2 V2 Xu 2 V 2 ~ UU 2 V2 

U 2 = V 2 "~ r~7J — , V 2 = rz ; , 

Xv 2 - UV 2 U 2 Xv 2 - \1V2 
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the form of which coincides with the lattice modified Boussinesq system ([79]). 



Theorem 5 The linear system 



d a ^ = K n (u, ui,9 a u; a)* 



(34a) 



where 



( 



(V 



(nu—adau) a\(avd a u+u(nv—ad a v)) 



K n (u,ux,d a u;a) : 



: 3 -A 3 



1 



V 




(34b) 



constitutes a Lax pair of system B nm , In other words, the compatibility condition 

dpK n {a, ui,9 a u;a) - d a K m (u, u 2 , dpu; 0) + [K n (u, ui, <9 Q u; a), K m (u, u 2 , dpu; (3)} = 
holds on solutions of system (132\) . 

5 Conclusions 

In this paper we have investigated the symmetry structure and conservation laws for two key examples of higher 
order integrable discrete systems, namely the lattice Boussinesq and the lattice modified Boussinesq system of 
partial difference equations. A third example in this class is the lattice Schwarzian Boussinesq equation, [9], which 
can also be written as a quadrilateral system of three fields, ifTTl . and for which we expect the symmetry analysis 
to be similar to that of the one performed in Section [2] Last year, in J4l, a systematic search for integrable cases 
of 2- and 3-field systems of Boussinesq-type form revealed the existence of certain parameter-deformations of the 
Boussinesq lattice systems, and recently in |[2H it was shown that these systems arise from a direct linearization 
scheme similar to the one exploited in Ifl2l but involving a more general dispersion structure. We expect that 
the deformed Boussinesq systems, including the analogues of the Schwarzian Boussinesq case, admit a similar 
analysis as the one performed in the present paper, but we postpone those generalizations to a future publication. 

Another related aspect, which is currently in preparation, is the existence of Lagrange structures for the Boussi- 
nesq systems. The scalar 9-point equation £[]) admits a natural Lagrangian, which was already reported in Ifl2l . but 
a Lagrangian for the modified lattice Boussinesq equation © was not given. In fact, in Q it was shown that the 
entire lattice Gel'fand-Dikii hierarchy, which comprises the lattice potential KdV equation as well as the lattice 
Boussinesq system and associated higher order systems, admits a Lagrangian structure which can be written in a 
surprisingly compact form. Recently, Lagrangians were established for the modified Boussinesq system compris- 
ing the fully discrete, semi-discrete systems as well as the corresponding generating PDE, |6]. These Lagrangian 
structures allow us to investigate the symmetries and conservation laws for these integrable systems from the point 
of view of the relevant Noether theorems. These aspects are part of ongoing research. 
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